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We consider the four-boson and 3+1 fermionic problems with a model Hamiltonian which encap- 
sulates the mechanism of the Feshbach resonance involving the coherent coupling of two atoms in 
the open channel and a molecule in the closed channel. The model includes also the pair-wise direct 
£S) , interaction between atoms in the open channel and in the bosonic case, the direct molecule-molecule 

interaction in the closed channel. Interactions are modeled by separable potentials which makes it 
possible to reduce the four-body problem to the study of a single integral equation. We take ad- 
vantage of the rotational symmetry and parity invariance of the Hamiltonian to reduce the general 
eigenvalue equation in each angular momentum sector to an integral equation for functions of three 
real variables only. A first application of this formalism in the zero-range limit is given elsewhere 
[Y. Castin, C. Mora, L. Pricoupenko, Phys. Rev. Lett. 105, 223201 (2010)]. 



^ t Resume : Equations integrates pour le probleme a quatre corps 



Nous considerons les problemes a quatre bosons et a 3+1 fermions en utilisant un Hamiltonien 
modele incluant le mecanisme de la resonance de Feshbach avec un couplage coherent entre deux 
atomes dans la voie ouverte et une molecule dans la voie fermee. Le modele comprend aussi une 
interaction directe entre atomes dans la voie ouverte, et dans le cas bosonique, l'interaction entre 
deux molecules de la voie fermee. Les interactions sont modelisees par des potentiels separables, ce 
qui permet de ramener le probleme a quatre corps a la resolution d'une equation integrale unique. 
L'invariance par rotation et par parite du Hamiltonien permet de reduire l'equation aux valeurs 
propres generale dans chaque secteur de moment cinetique fixe a une equation integrale portant sur 
des fonctions de trois variables reelles seulement. Une premiere application de ces considerations a 
deja ete mise en ceuvre dans [Y. Castin, C. Mora, L. Pricoupenko, Phys. Rev. Lett. 105, 223201 
(2010)]. 
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O . I. INTRODUCTION 

The few-body problem is the object of a renewed interest considering the possibility to study experimentally 
this problem in a resonant regime with cold atoms close to a Feshbach resonance 0. In this resonant regime, the 
s-wave scattering length a describing the interaction between the particles is much larger (in absolute value) than the 
range b of the interaction potential, which may be estimated by the van der Waals length of a few nanometers [1| . 

The three-body resonant problem, initiated in particular in Refs. jil-Q, bas been the subject of numerous studies 
[l|. Two distinct situations have been identified, depending on the quantum statistics and the mass ratios of the 
particles. In the first case, there exists a zero-range limit, where both the true range b and the effective range r e of 
the interaction tend to zero. The interactions may then be replaced by contact conditions on the wavefunction, which 
define the so-called Bethe-Peierls @, Q or zero-range model. The interaction is then described by a as the single 
• • , parameter. The paradigm of this first case is three spin- 1/2 fermions with equal masses for the two spin states. In 
the second case, the EEmovian case, the zero-range limit does not exist and is replaced by a limit cycle. The minimal 
model is then a modified zero-range model supplemented with three-body contact conditions involving the three-body 
parameter R t that depends on the microscopic details of the interaction. Then the interaction is described by only 
two parameters: a and Rt. The archetype of the second case is three indistinguishable bosons. The second case leads 
to the Efimov effect, that is for \a\ — oo there is an infinite number of weakly bound trimers with an accumulation 
point in the spectrum at zero energy. 

In reality, the physics of the Feshbach resonance is not exhausted by the zero-range model In particular, a 
Feshbach resonance involves two channels in the interaction between two low-energy atoms, one channel being open 
and the other one being closed due to energy conservation. In the specific case of ultra-narrow Feshbach resonances, 
where the effective range r e is negative and much larger in absolute value than the true range b, r e becomes an 
important parameter. It may be taken into account by modified Bethe-Peierls contact conditions @4ll|j the so-called 
effective range approach, or more straightforwardly by use of a two-channel model [T2-20]. In the general case, 
Feshbach resonances are richer, as e.g. a form or potential resonance may take place in the open channel, in which 
case the background scattering length o,b g is large and the modified Bethe-Peierls conditions [§4lll may not reproduce 
correctly the momentum dependence of the two-body scattering amplitude at low energy [41| [20l - |23j . Using a two- 
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channel model, with the description of the direct atomic interaction in the open channel by a separable potential, is 
then the simplest way to have a reasonably accurate description of the atomic physics behind the Feshbach resonance, 
while keeping the simplicity of the integral equations resulting from zero-range models. 

The four-body problem is now under active investigation. The situation is more open, since no analytical solution is 
available (contrary to the three-body case). In the absence of Efimov effect, the plain Bethe-Peierls model may be used 
for broad Feshbach resonances [IS, HE) • In principle exact numerical calculation of the resulting integral equation may 
be done, and was indeed performed to determine a dimer-dimer scattering length in Ref. [25| . Also efficient numerical 
methods have been used to solve the case of four or more harmonically trapped two-component fermions [26W28I ] In the 
Efimovian case, for four bosons, finite range real space models, specified by an interaction potential V(rij), have been 
recently used in the context of cold atoms but have been solved only in an approximate way 29] . As a consequence, the 
tetramer states predicted in (29| that have (negative) energies larger than the trimer ground state energy are actually 
resonances rather than genuine stationary states. In the nuclear physics context, however, exact numerical solutions 
of the four-body problem with real space interaction potentials V(rij) can be obtained thanks to the reformulation 
in terms of momentum-space Faddeev-Yakubovsky integral equatio ns pOl l3lj | . A numerical solution of a low-energy 
effective theory based on this integral formulation was used in Ref. [32j to predict universal properties of the ground 
four-boson state of that theory. The Faddeev-Yakubovsky integral formulation may also be used to solve the four-body 
problem with a separable potential interaction and to calculate the energy width of the resonances corresponding to 
the tetramers of [20|, as recently done in [33|. A limitation of the present Faddeev-Yakubovsky formulation seems to 
be, however, that it does not encapsulate the atomic physics of the Feshbach resonance. 

To answer as quantitatively as possible current questions on the four-body problem, we use in this paper a four- 
parameter two-channel model to derive without approximation a simple integral equation in momentum space for 
the four-body problem, involving an unknown function of two atomic momenta only. No solution of the equation is 
presented here, but we presented elsewhere a first application of this formalism in the limiting case of a zero-range 
interaction [34| and we expect that other applications will come. In this paper, we take as specific examples first 
the case of four spinless bosons, in section [H] and then the so-called 3 + 1 fermionic problem of three same spin 
state fermions of mass m interacting with a distinguishable particle of mass M, in section IIIII The reduction of 
the integral equation thanks to the rotational symmetry of the Hamiltonian is performed in section IIV1 leading to 
unknown functions of three real variables only. We conclude in section [V] 

II. THE FOUR-BOSON PROBLEM 
A. Model Hamiltonian and two-body scattering amplitude 

Our model is a two-channel model of a magnetic Feshbach resonance [l2|, [l3[ . This means that the atoms may 
actually exist in two different forms in the model, either as atoms in the open channel, treated as elementary bosons, 
or as molecules in the closed channel, the so-called closed channel molecules. Since these molecules have a size of the 
order of the van der Waals length b, usually smaller than the mean distance between them, they are also treated as 
elementary bosons. There exists a short-range coupling between the two channels, due in particular to the existence 
of hyperfinc atomic degrees of freedom, so that two atoms may coherently be converted to a closed-channel molecule, 
and vice versa. This coherent interconversion leads to an effective interaction between the atoms, that may be made 
resonant by tuning with a magnetic field the energy E mo \ of a closed-channel molecule counted with respect to the 
dissociation limit of the open channel potential. Even in the absence of this interconversion, atoms in the open channel 
can interact via the van der Waals potential, and this so-called background interaction is sometimes also resonant 
and cannot be neglected. 

Our model Hamiltonian thus involves two coupled bosonic fields, one for the open channel atoms and one for the 
closed-channel molecules. For simplicity, we assume that there is no trapping potential and we perform calculations 
directly in momentum space without any quantization volume (the real space version of the Hamiltonian was written 
in Ref. [35|). The creation 4 and annihilation ak operators of one atom of wave vector k, and the creation 4 and 
annihilation &k operators of one closed channel molecule of ccntcr-of-mass wave vector k then obey the canonical 
commutation relations 

[Ok, 4'] = [ & k, 4] = (27T) 3 £(k - k'). (1) 

The other commutators are zero, in particular ak and 4 commute with &k and b k . The Hamiltonian may be split in 
four terms, 



H — Hat + Hmol + H opcn + H at 

— mol- 



(2) 
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The first two terms contain the kinetic energy (and the internal energy) of non-interacting atoms and molecules, 

d 3 k 



rE^a^a^ 



(2tt) 3 ' 

H mo \ = J j^—j^{E mo i + -Ek)bl.bk. 
We have introduced the convenient notation for the kinetic energy of an atom of wave vector k, 

h 2 k 2 



Eu = 



2m 



(3) 
(4) 

(5) 



and we shall use in what follows the parity invariance E^ = .E-k- 

The third term in Eq. ([2]) describes the direct interaction between atoms in the open channel. In princip le, one 



should describe this interaction using true binary potentials Vfcj) with a 1/Vy Van der Waals tail [36H38I ]. This 
however precludes the derivation of a simple inte gra l equation for the four-body problem. We thus rather take the 
convenient form of a separable potential, see e.g. [391 ]: 
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d 3 ki d 3 k 2 d 3 k 3 rf 3 /c 4 k 3 — k 4 k 2 
(2tt) 3 (2tt) 3 (2^)3 (2^)3 2 



ki 



)(2TryS(k 1 + k 2 - k 3 - k 4 )a^ a k a k2 a kl 



(6) 



It involves a bare coupling constant go and a cut-off function \ making the theory well defined. In practice, we take 
as in Ref. [H], EH a real Gaussian cut-off function, which is both a smooth and rapidly decreasing function of k, 



X(k) 



-k 2 b 2 /2 



(7) 



where b is of the order of the van der Waals length. The Gaussian choice is of course arbitrary, and may be replaced 
by any rapidly decreasing function of k. An important point is that the separable potential is able to capture the 
essential features of the Feshbach resonance, such as the dependance of the scattering length with the magnetic field 
(fT2")l . or the structure (fT5|) of the effective range. 

The last term in Eq. @ describes the coherent interconversion of pairs of atoms into closed-channel molecules, 
responsible for the Feshbach resonance: 



H, 



at — mol 



= A 



d 3 h d 3 k 2 
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kj 



) & k 1 +k 2 «k 2 ak 1 +x(- 
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)« kl « 
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It involves a coupling amplitude A between the two channels, taken to be real, and (for simplicity) the same cut-off 
function \ of the relative momentum of two bosons, implementing in momentum space the fact that, in real space, the 
conversion process is non-local over a radius sa b only, that is two atoms are converted to a closed-channel molecule 
with an appreciable probability amplitude only if their interatomic distance is of the order of b or less. This cut- 
off function is actually required to make H at -mo\ mathematically sound. Note that each of the four terms of the 
Hamiltonian conserves the total momentum and preserves Galilean invariance, as it should be. 

It appears that our model Hamiltonian depends on four bare parameters b, go, E mo \ and A. Apart from the true 
range b, the other three parameters are more physically expressed in terms of the scattering length a, the background 
scattering length db g and the always positive Feshbach length [l(| ■ The link is performed thanks to the solution of 
the two-atom scattering problem, that was detailed i n 13 51 for two opposite spin fermions, and that is adapted to the 
present bosonic case by replacing A 2 in the result of [35J by 2A 2 . The resulting expression of the two-body T matrix 
is then 



Airfr 2 

k f \T(E + l 0+)|k j; ) = X(k/)x(k i ) f(E + z0+) 



(9) 



where k^ and k/ are the wave vectors of the incoming wave and outgoing wave in the center of mass frame, E is the 
energy in that frame and we have defined the function of energy f(E + i0 + ) such that: 



4irh 2 f(E + i0+) 



2A 2 
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E - E n 



d 3 k 



X 2 (k) 



(2tt) 3 E + i0+ - 2E k 



(10) 



The calculation in [35J was done on shell, that is for k 2 = k 2 = mE/ti 2 , and E > 0. In this case, the function f(E) is 
simply related to the s-wave scattering amplitude fk at that energy, with k — (mE) 1 / 2 /h, by the relation 

/fc = x(k/)x(k 4 )/(£ + *0+). (11) 
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By an analytic continuation argument, Eq. ([9]) actually also holds off shell, for all values of k^, k/ and E. Assuming 
that E mo \ is an affine function of the magnetic field B, a good approximation if B remains close to the Feshbach 
resonance location Bo, one extracts from Eq. ([9]) the standard formula for the magnetic field dependence of the 
scattering length a = limfc^o — 



1 - 



AB 

B — Bo 



(12) 



where db g is the background scattering length and AB is the width of the Feshbach resonance. The Feshbach length 
is connected to this width by 



R,. 



ibgSfiAB 1 



(13) 



where 6/j, = dE mo \/dB. The physical parameters of the Feshbach resonance are finally found to be related to the bare 
ones by 
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JL = 



1 



47rS 2 /m 

9o by/n 

Airh 2 /m 
9o - 2A 2 /E mol ~* 

9l 



b^ 



8ttA2< 



(14) 
(15) 
(16) 



which allows a direct connection with the experimental parameters. The expression of the function f(E) of our model 
for E > 0, and of its analytic continuation to E < 0, is also quite simple in terms of the physical parameters (35| . 



f(E + i0+) 
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fk 



^b A 



1-1- 



flbti 



k 2 



k 2 -Q 2 



ikerf (—ikb) + ik 



(17) 



where erf is the error function, we have set k = (mE) 1 / 2 /h for E > and k — i(-mE) 1 / 2 /h for E < 0, and we have 
introduced Q 2 such that — 1/Q 2 = ab g -R*(l — a^g/a). This allows in particular a direct calculation of the effective 
range r e of the model (35|, defined as —l/fk = -+ik— \k 2 r e + o(k 2 ) for k 0, 



Wi-^) 2 + it-^. 

V a / \/tt 



262 
a 



(18) 



Another quantity of interest for what follows is the internal energy E mo \ of a decoupled closed-channel molecule. 
It may be expressed in terms of b and of the three physical parameters a, Ob g and i?* using Eqs. (|14I15I16|) . For 
simplicity, we give this expression only at resonance 1/a = : 



TTTCS 



mi?* (Obg — by/it) 

Since i?„ is positive, -E^oi ^ s negative if and only if ab g > b\pn. 



(19) 



B. Four-body ansatz and coupled equations 



We consider for simplicity a four-body eigenvector of our model Hamiltonian in free space with an eigenenergy E 
obeying the two conditions 

E < and E — 2E mol < 0. (20) 

The case of E > is a four-atom scattering problem that may be treated along the lines of e.g. [l8j], including in 
particular an incoming free four-atom state. For the same reason of simplicity, to avoid the treatment of the scattering 
problem of two closed-channel molecules, we also impose E — 2E mo i < [42[ In the two-channel model, the four-body 
state vector is in general a coherent superposition of a component with four atoms in the open channel, a component 
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with two atoms in the open channel and one closed-channel molecule, and a component with two closed-channel 
molecules, 



I*) = 1^4 at) + IV^at+lmol) + \tp2 mol)- 

Restricting without loss of generality to a zero total momentum state, we obtain the ansatz 



1^4 at) = 
|^2at+lmol) = 



2 mol / 



d 3 k x d 3 k 2 d 3 k 3 d 3 k 4 
(2tt) 3 (2tt) 3 (2tt) 3 (2tt) 3 
d 3 fci d 3 k 2 t 

d 3 k , , 

(^( k )^ T - k |o> 



(2^) 3 ,5(k 1 + k 2 + k 3 + k 4 )A(k 1; k 2 , k 3 , k 4 )4 1 oJ ca aJ C3 oJ C4 |0) 



,t „t 



(k!+k 2 )"k 
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(21) 

(22) 
(23) 
(24) 



Thanks to the bosonic symmetry, ensuring that all the creation operators commute, we can impose that the amplitude 
A is a symmetric function of its four arguments, the amplitude B is a symmetric function of its two arguments, and 
C(k) is an even function of k. 

It remains to calculate H \^) and to project Schrodinger's equation = (H — E)\^f) on the three relevant subspaces. 
Projection in the sector with four open channel atoms yields 







d 3 k± d 3 k 2 d 3 k 3 d 3 k 4 
(2tt) 3 (2tt) 3 (2tt) 3 (2tt) 3 



(27r) 3 <5(ki + k 2 + k 3 + k 4 )F 4 (ki, k 2 , k 3 , k 4 )4 1 4 2 al C3 4jO) 



with the function 

i^(ki,k 2 ,k 3 ,k4 



-E + J2 E K, 



A(k 1 ,k 2 ,k 3 ,k 4 ) + x(^ L ^ i )[ 6 ffoS(k 3 ,k 4 )+Ai?(k3,k 4 )]. 



(25) 



(26) 



The first term in the right-hand side clearly contains the kinetic energies of the four atoms. In the second term, 
the occurrence of the B amplitude of the ansatz is due to the conversion of a closed-channel molecule into two open 
channel atoms, as is revealed by the presence of the factor A. Finally, the term proportional to 50 is due to the direct 
pair-wise interaction between the four atoms in the open channel. Since this interaction is modeled by a separable 
potential, it involves a function of two momenta only, obtained by contracting the four-boson amplitude A with the 
cut-off function x, 



S(ki,k a ) 



' /!/ '' : ' / ' ;fc4 -(2^(k 1 +k 2 +k 3 + k 4 ) x (^-^)A(k 3 ,k 4 ,k 1 ,k 2 ). 



(2tt) 3 (2tt) s 



(27) 



We have introduced the notation B to underline the formal similarity of this term in Eq. (126)) with the two-channel 
contribution KB. Eqs. (|25|26p allow to express the four-boson amplitude A in terms of two-body amplitudes B and 
B. Because of bosonic symmetry, Eq. ([23]) does not imply that F 4 is identically zero, it only implies that the totally 
symmetric component of F 4 , symmetrized over the 4! permutations er of the four arguments ki !2j3j4 , is zero. This 
leads to 



t t vl 1 \ - ^, k CT(1) -k j(2) 6.90^( ^(3), k CT(4) )+Ag(k CT(3) ,k CT(4) ) 
A(ki,k 2 ,k 3 , k 4 j = — X( 7. ) 



o-es 4 



(28) 



The unknown amplitude A can thus be eliminated in the remaining part of the calculation. Insertion of Eq. 
the definition (|2"T|) gives 



S(ki,k 2 ) = 



d 3 k 3 d 3 k 4 (2i:) 3 S(j: 4 n =i k ») 1 
(27r) 3 (2^) 3 E-Y: 4 n=1 E kn 4! 



£*(- 

cr£S 4 



-)x( 



K(i 



(i) 



C CT (2) 



)[6goB + AB}(K {3) ,K {4) ). (29) 



Collecting permutations giving identical contributions, up to a relabeling of the integration variables, gives the more 
explicit form [43[ 



B(ki,k 2 ) = 



d 3 k 3 d 3 k 4 (27r) 3 (5(^ =1 k n ) 1 ,k 3 -k 



+x( 



(2tt) 3 (2tt) 3 E-J2l =1 E k (> 
ki - k 2 



«x(- 



-)[6g B + AB](k u k 2 ) 



-)[6. 9o S + AB](k 3 ,k 4 )+2 



) x(- 



x( k3_ki )[65o ^ + AB](k2;k4) + 1 2 



(30) 



6 



This equation is the first important equation of this subsection. 

We now project Schrodinger's equation = (H — E)\*5>) in the sector with two open channel atoms and one 
closed-channel molecule. We obtain 

with 

F 2 (k l7 k 2 ) = [ J B kl + J B k2 + i£ kl+k2 + J B mol - J B] J B(k 1 ,k 2 ) + 12A J B(k 1 ,k 2 ) 

+ -)[2AC(ki + k 2 ) + gopfa + ka)]. (32) 

In the first term in the right-hand side of F2, we recognize the kinetic energies of two atoms, plus the kinetic and 
internal energies of a closed-channel molecule of wavevector — (ki + k 2 ). The contributions KB and AC are due to the 
interchannel coupling, respectively converting a pair of atoms into a closed-channel molecule and vice versa. The fact 
that the same contraction B defined in Eq. (|2"T1) appears as in the four-atom sector is due to the choice of the same 
cut-off function in the separable potential and the open-to-closed-channel coupling. Finally, the term proportional 
to go i n Eq. (|32[) results from the direct open channel interaction between the two open channel atoms present in 
|^2at+imoi)- This involves the contraction of B with the cut-off function x, that we call /3, 

0(K) = J ^^(2^(k 1 +k 2 -K) x ( k l T ^)B(k 1 ,k 2 ). (33) 

Since F2 is a symmetric function of ki and k 2 , Eq. (|3Tj) directly imposes F 2 (ki, k 2 ) =0. 

Finally, projecting Schrodinger's equation = (H — E)\ty) in the sector with zero open channel atom and two 
closed-channel molecules gives 

= / (^3^o(k)6 k 6L k |0) (34) 

with 



F (k) = (2E mol + £ k - E)C(k) + A[/3(k) + /3(-k)] 



(35) 



The first term in the right-hand side contains the kinetic and internal energies of the two closed-channel molecules, 
and the last term originates from the conversion of two atoms into a closed-channel molecule. The direct open channel 
interaction does not enter here. The fact that the same contracted function (3 appears as in Eq. (|3"2"f is due to the 
fact that the same cut-off function \ is used in H opcn and -ff a t-moi- Since -Fb(k) is an even function of k, Eq. (|34[) 
directly imposes Fo(k) = 0, which allows to eliminate C, 



C(k) = A 



/3(k)+/3(-k) 
E — Ek — 2E mo \ 



(36) 



We used the condition Eq. (I20[) to ensure that the energy denominator in that expression cannot vanish. Reporting 
this into the equation F2 = gives the second important equation of this subsection, 



12AS(ki,k 2 ) = [E kl +Eu 2 + -E Wl+ u 2 +E mo i-E]B(k u k2) 



+ x( 



ki 



2A 



2 /3(k 1 +k 2 )+/3(-k 1 -k 2 



E - -E kl + k2 



2E, 



mol 



(37) 



C. A single integral equation 



We now show that, remarkably, the rather involved set of equations (f30|) and (137)) can be reduced to a single explicit 
integral equation on the unknown function _D(ki,k 2 ) defined as: 



A£>(ki, k 2 ) = 6 3o 5(ki, k 2 ) + A£(ki, k 2 ). 



(38) 
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We note that this function naturally appears in the right-hand side of Eq. (|30|) . The left-hand side of Eq. (|30[) involves 
B, which is a linear combination of the functions D and B. The only non trivial step is thus to express B as a function 
of D using Eq. ff37|) . We rewrite Eq. ([37]) introducing an effective (3 function, 



ft.(K )sW K) + ^ K ' + «- £ K > , , 39 ) 
and introducing what we call the relative energy for a reason that will become clear in the final result, 

£ ro i(ki, k 2 ) =E-E kl - E k2 - ^ kl+k2 . (40) 
Eliminating B in terms of D and -B thanks to Eq. ([38]) . we transform Eq. (l37l) into 



£>(ki,k 2 )_ Dn . . ff 2 .ki-k 2 ./3 off (ki+k 



u(ki,k 2 ) v " z; 2A 2AV 2 ' u(ki,k a ) 

with the notation 

.9o 



(kx.ka) = 1 - ^ [S mo] - ^(kx.ka)]. (42) 



Eq. (|4T|) docs not immediately give £> in terms of D because of the occurrence of (3 c s in the right-hand side. Conse- 
quently, one has first to express /3 e ff in terms of D. We note that the Hamiltonian H of the system is invariant by 
parity. We can thus without loss of generality assume that the eigenstate | of H that we consider is of well defined 
parity 77, with 77 = 1 (even state) or —1 (odd state). Hence, the functions B, B and thus D, f3 have also the parity 77. 
For example, 

D(-ki,-k a )=T/I>(ki,k 2 ) and /?(— K) = r//3(K). (43) 

This makes Eq. (l39|) local in K, expressing /3 e g(K) as the product of /3(K) and of a known 77-dependent function 
of K. Since /3(K) results from the contraction of _B(ki,k 2 ) with the cut-off function x, see Eq. (1331) . it remains to 
multiply Eq. (gTJl by <5(ki + k 2 - K)x ( kl ~ k2 ) and to integrate over k x and k 2 . In the left-hand side of Eq. (gTJ| the 
contraction of D/u with x immediately appears, we call it 

/" d 3 k x d 3 k 2 , , 3 , . ki - k 2 L>(ki,k 2 ) 

7 (27r)3 (27r) J 2 u(ki,k 2 ) 

We finally obtain the desired expression /3 e g in terms of D, 

(3 cS (K) = «(K)A(K) (45) 

where we have introduced the parity dependent notation 

(1 + 77)2A 2 /.g T 1 ff0 2 /• d 3 fc X 2 (k) 1"' 



v(K) 



1 + 



E-Ek- 2E n 



2A 2 J (2tt) 3 u (|K - k, |K + k) 



(46) 



Eq. (|4ip thus now gives S explicitly as a functional of D. This, combined with Eq. ([58)1 . immediately gives B as a 
functional of £). The left-hand side of Eq. (|30l) is then also expressed in terms of D. Since the right-hand side of 
Eq. (|30p is directly a functional of D, Eq. (|30| (multiplied for convenience by the factor 6/ A) reduces to the desired 
closed integral equation for the unknown function D. Collecting in the right-hand side of that integral equation the 
two terms involving -D(ki, k 2 ), that is one term originating from B and the other term originating from the right-hand 
side of Eq. ([3"0")l. we get as a prefactor of Z)(ki,k 2 ) the function 

Praetor = -1 ( I - -J—) + / ^ ^(^1). (47) 

go \ u(k u k 2 )J J (2^)3(2^)3 E-Yt =l E* n 

This may be transformed using the explicit expression of the function u, the identity 

(^) 3 S(ELi k ») _ (2^) 3 <5(ki + k 2 + k 3 + k 4 ) 

E - En=l ^k„ Srel(kl,k 2 ) - 2^ (k4 _ k3)/2 



(48) 
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and performing the change of variable of unit Jacobian, k3. 4 = ^ ± k. Remarkably, the prefactor can then be linked 
to the function / appearing in the two-body T matrix at the energy E ie \, 



prefactor 



47T/l 2 /[£rcl(kl,k 2 )]- 



(49) 



This justifies the name relative energy for E TC \. To conclude, we give the final equation for D, which is the main result 
of this section: 



= 



mD (ki,k 2 ) 
4 7 r? l 2 /[£; rcl (k 1 ,k 2 ) ' 2A 

k 3 - k 4 



go ki -k 2 ,«(ki +k 2 ) 
TTTa *( = ) TIT. ,. \ A ( k i + k 2) 



u(ki,k 2 ) 



d 3 k 3 d 3 h J (2nfS(Y:Li kn) 



J (2tt) 3 (2tt) 3 I 
k 2 



x( k ^) J D(k 3) k 4 ) + 2 X ( k ^i) J D(k 2 , k 4 ) + 2 X (^ 



)£>(k 1) k 4 ) 



(50) 



This has to be solved for a bosonic symmetry of D, £>(ki,k 2 ) = Z?(k 2 ,ki). We note that the second term in the 
right-hand side of Eq. ([50)1 has no equivalent in zero-range theories and is to our knowledge new. 



D. Including closed-channel molecule interaction 

In the previous subsections, any direct interaction between closed-channel molecules was neglected. Here we derive 
a modified closed integral equation for the function _D(ki,k 2 ), including the intermolecular s-wave interaction in the 
form of a separable potential with the previously introduced cut-off function x(k) and a bare coupling constant g™ . 
This amounts to adding to the model Hamiltonian Eq. ^ the term 

H 2mo i = \gT X J ^^^3^^x(k)x(k0fo t i K+k ,& t i K _ k ,6iK+k^K-k- (51) 



This term i? 2mo i has a non-zero action only on the component iV^moi) of the four-body state vector \^f) in Eq. pip . 
Projecting Schrodingcr's equation in that two closed-channel molecule sector gives Eq. (|34| with the function i"b(K) 
modified as 

F (K) = (2E mol + E K - E)C(K) + A [0(K) + /3(-K)] + g^° l J 0^ x (K) x (K')C(K'). (52) 

Since Fo is an even function, Schrodinger's equation requires Fq(K) = 0. This leads to an integral equation of the 
Lippmann-Schwinger type: 

C(K) = C (K) + 5o mol G° 1 oi(K)x(K) J — X (K')C(K'), (53) 

where C (K) = AG° nol (K) [(3(K) + /3(-K)] , and G^ ol (K) = 1/(J^ - E k ) is the free Green's function of the two 
closed-channel molecules at the relative energy E™^ = E — 2E mo \. We recall that here E < 2E mo i, see Eq. (|20|) . To 
solve the Lippmann-Schwinger equation (|53[) . one multiplies it by x(K) and one integrates over K. A closed equation 

is then obtained for 7 = J {* 2 J^ 3 x(K)C(K). This leads to the expression of C(K) as a functional of Cb(K), and thus 

of the contracted pair function /3(K): 

C(K) = Co(K) + g mol G? nol (K)x(K) J X (K')C (K'), (54) 

where g m ° l is the effective molecular coupling constant defined by 

1 t d 3 K' 

K gf^~J w? 

Similarly to the previous subsection, the four-body problem at an energy satisfying Eq. (|20[) can be reduced to 
an integral equation for the function Z)(k 1 ,k 2 ) defined in Eq. (|38[) . Remarkably, one finds that the only change to 
perform to Eq. (|50[) is the substitution 

2A 2 r d 3 K' x(K')G^ ol (K')A(K') 



9 nWl - ( 3- / ^G? nol (K'k 2 (K') ) • (55) 



( 2A f d K 

1 ff mol x(K>(K) J ~lfJ {2*r n.(K')i(K>) y 
v(K)A(K) -> u(K)A(K) + — — <{ - 2(1+ „ )A2gmol rf3g , [x (kq GL( kt > ( 56 ) 

J R(K' 



1 + 
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where we have introduced the functions: 



i + 2(1 + ^?)A 2 Q . r) 2A 2 _ f d 3 k 



X 2 (k) 



(57) 



As expected, the odd states (77 = —1) are not affected by the s-wave closed-channel molecule interaction. For even 
states (77 = 1), the term added by the substitution is in general non zero. It may significantly affect the integral 
equation for D, that is even for the low-energy eigenstates \E\ <C H 2 /(mb 2 ), when the following resonance condition 
is met: 



1 



4A 2 g 



2„mol 



fit) 



d 3 K' [ X (K')G° mol (K')} 



(2tt) 3 



R(K') 



< 1. 



(58) 



To obtain Eq. (|55|) . we had to recalculate the function /3 e g = /3(K) + (2A/ g )C(K), now including the closed-channel 
molecule interaction. For this purpose, we used the expression of A in terms of /3 and /3 c g: 



A(K)=/3(K)-/MK)x|,/||^ 



X 2 (k) 



K 



k,±K 



k) 



Introducing as in the previous subsection the parity 77 of the state one can write: 



MK) 



2(l + r7)A 2 G° nol (K 
9o 



/3(K)+ 2(1 + ??)A2gm0l G° nol (K) x (K) / ^ x(K ')G? nol (K')/3(K') 
5o 7 (2tt) 3 



(59) 



(60) 



Simple algebra on this expression and comparison of the result to Eq. (|4"5"j) finally lead to the prescribed substitution 
Eq. (151)1) . An interesting extension of this subsection would be to include a direct interaction between the closed- 
channel molecules and the open-channel atoms, to see if there still exists a closed integral equation for D(ki, k^). 



III. THE 3+1 FERMIONIC PROBLEM 



We now consider the four-body problem corresponding to the interaction of three same spin state fermions of mass 
m with an extra, distinguishable particle of mass M, for example of another atomic species. There is no possible 
s-wave interaction among the fermions, and we assume that their interaction in odd angular momentum waves (in 
particular the p-wave) is negligible. The fermions thus interact only with the extra particle, either directly in the open 
channel, or indirectly through the creation of a closed-channel molecule. Since the extra particle is alone, there is no 
need to specify the statistical nature of that particle, nor of the closed-channel molecule. An important difference with 
the previous four-body bosonic case is that the subspace with two closed-channel molecules will not be populated, and 
also that there is no direct open channel interaction in the subspace with two atoms and one closed-channel molecule. 
All this simplifies the problem. 



A. Model Hamiltonian and two-body scattering amplitude 



The fermionic annihilation and creation operators are called Ck and c^, they obey the usual free space anticommu- 
tation relations {ck, c^,} = (27r) 3 5(k — k'). They are assumed to commute with the extra particle creation operator 
and with a closed-channel molecule creation operator b^, that also commute among themselves. The total Hamiltonian 

H = -Hat + -Hmol + ffat-mo] + H opcn with 



H mo \ 

Hg^t— mol 
-Hon on 



.90 



d 3 k 
(2tt) 3 

d 3 k 
(2tt) 3 



E mo i 



l + a 



E k } & k 6 k 



A / d [(frr) 3 P *( kl2 )[ & ki+k 2 a kiCk 2 + « kl c kA 1+ k 2 ] 
d 3 k 1 d 3 k 2 d 3 k 3 d 3 k 4 



[(2tt; 



:m 



X(ki 2 )x(k43)(27r) 3 (5(ki + k 2 - k 3 - k 4 )aj C4 4 3 Ck 2 a kl 



(61) 
(62) 
(63) 
(64) 
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We have introduced the free fermion dispersion relation E^ = h 2 k 2 /(2m) and the mass ratio of a fermion to the extra 
particle: 

" s w < 65 > 

As compared to the bosonic case, we have dropped the factor 1/2 in front of go in H opcn since the extra particle is 
distinguishable from the fcrmions. More important, to maintain the Galilean invariance, the argument of the cut-off 
function x is now the relative wave vector of a fermion (of momentum hk 2 ) with respect to the extra particle (of 
momentum Kki): 

/k 2 ki\ k 2 -aki 

ki 2 = /i — = — — , 66) 

\m M ) 1 + a 

where /i = mM/(m + M) is the reduced mass. In this notation for the relative wavevector ki 2 , it is important to 
keep in mind that the first index refers to the extra particle and the second index to a fermion. 

The two-body scattering of the extra particle with a fermion is similar to the two-body bosonic problem of subsection 
III Al In the bosonic formulas Eqs. (|14ll5ll6p . one simply has to replace m with 2/i and 2A 2 with A 2 . The two-body 
T matrix is now given by 

(k f \T(E + i0+)|k i ) = -^! x ( k/ ) x (k 4 ) f(E + *0+) (67) 

A* 



where k^ and k/ are relative wavevectors and the function / now reads 



2nh 2 f(E + 



A 2 

go + E—R 



mol 



d 3 k x 2 (k) , , 

v (68) 



(2n) 3 E+10+ -h 2 k 2 /(2^)' 



The connection of / with the s-wave scattering amplitude ff. on shell is still given by Eq. (ITT|) with now k = fc, = 
k f = {2^E) 1 ' 2 /h. Also relation Eq. JX?) still holds, with k = (2^E)^ 2 /h for E > and k = i(~2fiE)^ 2 /h for E < 0. 



B. Four-body ansatz and coupled equations 

The ansatz for the four-body state is now simplified to |\&) = 1^4 at) + iV^at+i mol) with 



1^4 at) = 
1*02 at+1 mol) = 



rf 3 fc 4 d 3 k 2 d 3 k 3 d 3 k 4 
(2tt) 3 (2tt) 3 (2tt) 3 (2tt) 3 
d 3 h d 3 k 2 



(27t) 3 s(k 1 + k 2 + k 3 + k 4 )v4(ki, k 2 , k 3 , k 4 )4 1 4 a 4,4 4 1°) ( 69 ) 



(2tt) 3 (2tt) : 



-B(k u k 2 )bl 



(ki+k a )Mk 



Aio). 



(70) 



Taking advantage of the fcrmionic symmetry, we impose that A is an antisymmetric function of its last three vectorial 
arguments, and that B is an antisymmetric function of its two vectorial arguments. 

By projecting Schrodinger's equation = (H — E)\^f) in the four-atom sector, where we impose E < so as to 
avoid, as in the bosonic case, vanishing energy denominators, we obtain 







d 3 h d 3 k 2 d 3 k 3 d 3 k 4 
(2tt) 3 (2tt) 3 J2^f (2^ 



■(27r) 3 <j(k 1 + k 2 + k 3 + k 4 )F 4 (k l5 k 2 , k 3 , k 4 K i 4 2 4 3 4jo) 



(71) 



with 



^ 4 (ki,k 2 ,k 3 ,k 4 



-E + a£ kl + J2 E ^ 



A(ki, k 2 , k 3 , k 4 ) + x (k 12 )[3 5o -B(k 3 , k 4 ) + AB(k 3 , k 4 )]. (72) 



We have defined B by a contraction of A with the cut-off function x 

d 3 hd 3 k 2 



B(k 3 ,k 4 



[(2tt) 3 



X (k 12 )(2^) 3 <5(k 1 + k 2 + k 3 + k 4 )A(k l5 k 2 , k 3 , k 4 



(73) 
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Similarly to the bosonic case, the unknown function in the final integral equation will be a linear combination of B 
and of B: 

AD(k 3 , k 4 ) = 3g Q B(k 3 , k 4 ) + A£?(k 3 , k 4 ). (74) 

Since the antisymmetric part of F4 with respect to its last three vectorial arguments is zero by virtue of Eq. (|7ip , we 
obtain A as a function of D: 

A(k l7 k a , kg, k 4 ) = - A/3 — - [ x (k 12 )D(k 3 , k*) - x(k 13 )£>(k 2 , k 4 ) + x(k 14 )D(k 2 , k 3 )] . (75) 

E-(aE kl +2j„=2 £; k,J 

This we insert in the definition of B. After the renumbering 1234 — > 3412 in Eq. ([75)l . the extra particle now having 
the index 3, we also obtain B as a function of D: 



3~ f fksfk, (2n)36(ZLM , 1: v 



r)3]2 B-^+En^^)' 

x [ X (k 34 )£>(k 1) k a ) - x(k 31 )£>(k4, k 2 ) + x (k 32 ) J D(k 4 , k x )]. (76) 

By projecting = (H — E)\ty) in the subspace with two fermions and one closed-channel molecule, we obtain 

d 3 fci d 3 k 2 





with 



/ ^^ F ^^-^A<\°) (^) 



F 2 (k 1; k 2 ) = \ -E + E kx + E k2 + ^—E kl+k2 + E mol j J B(k 1 ,k 2 ) + 3AB(k 1 ,k 2 ). (78) 

Since F2 is an antisymmetric function, Schrodinger's equation directly imposes F% = 0. As in the bosonic case, we 
then introduce the relative energy 



£ rel (ki,k 2 ) = E — [ E kl + E k + -—E kl+k , , (79) 

1 + a 

and the function u(ki,k 2 ) = 1 — (go/A 2 )[£moi — ^rel(ki, k 2 )]. Then we obtain a second expression of B as a function 
of D: 

X^ k2 )=^(l-^y)^,k 2 ). (80) 



C. A single integral equation 

To obtain the desired integral equation on D, it remains to eliminate B in between Eq. ([76]) and Eq. (l80l) . Diagonal 
and non-diagonal contributions appear. In the diagonal contribution, to recover the two-body T matrix, we perform 
the usual change of variables of unit Jacobian, that introduces the relative wave vector k 34 and the center-of-mass 
wave vector K 34 = k 3 + k 4 of the extra particle of index 3 and one fermion of index 4: 

k 3 - — J— K 34 -k 34 (81) 
1 + a 

k 4 = ^_K 34 + k 34 . (82) 
1 + a 

Then, from the identity (1 + a)/m — X/fi, we obtain the simple relation 

(^) 3 5(E 4 n =i^n) (27r) 3 S(K 3 4 + k 1+ k 2 ) , , 

E-(aE k3 +j: n ^ 3 E k J E rel (k l! k 2 )-S 2 fci 4 /(2/x)' 1 ' 

Finally, the integral equation for D in the three fermions plus one extra particle problem is 
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fiD(ki,k 2 ) 



d 3 k 4 



(27T) 3 J(ELlkn) 



2irh 2 f[E rel (k l7 k 2 )] J (2tt) 3 (2tt) 3 £ - (a£ k3 + #k„) 



X (k 3 4)[x(k3i)^(k4,k 2 ) 



+ X (k 32 ) J D(k 1 ,k 4 )], (84) 



where the function f(E) is the fermionic one given by Eq. The integral equation (|84l) is the main result of this 

section. It has to be solved with a fermionic exchange symmetry for D, D(ki,k 2 ) = — Z)(k 2 ,ki). 

An instructive limiting case is to take the zero-range limit b — > in Eq. (|84p . assuming that a^ g and i?* are 0(b) 
whereas the scattering length a is fixed. Then the scattering amplitude of the model tends to the usual zero-range 
expression 



ik 



(85) 



In Eq. (|84[) one replaces the function \ with unity, and one performs the integral over the extra particle momentum 
k 3 . Setting E = — h 2 q 2 / (2/j,), q > 0, one obtains the integral equation for the zero-range Bethe-Peierls model for three 
same spin state fermions and one extra particle: 



1 + 2a 
(1~ 



2a 



(l + a) 2 

d 3 h 



2n 2 



-,1/2 



kik 2 



1 



->£>(ki,k 2 ) 

D(k 1 ,k 4 ) + J D(k 4 ,k 2 



k^ + kc^ + k 



la 
l+a 



(ki • k 2 + ki • k 4 + k 2 • k 4 



= 0. 



(86) 



This differs from the integral equation derived in [241 by numerical constants [44l | . The zero energy q — > and unitary 
1/a = limits of Eq. (|86| were recently used in [34| to study the emergence of a pure four-body Efimov effect (without 
three-body Efimov effect) in the system of three fermions plus one extra particle, as a function of the mass ratio a. 



IV. TAKING ADVANTAGE OF ROTATIONAL INVARIANCE 



In practice, integral equations such as Eqs. (|50I84[) are heavy to solve numerically because the unknown function 
a priori depends on 6 real variables. A significant reduction of the problem can be achieved by using the rotational 
symmetry of the Hamiltonian. This rotational symmetry implies the existence of degenerate eigenenergy subspaces of 
well defined total angular momentum I. Such a degenerate subspace is associated to an irreducible representation of the 
50(3) rotation group of spin I, and it is generated in Dirac's notation by 21 + 1 functions \l, mi), with \l, mi) of angular 
momentum mih along the arbitrary quantization axis z. The general solution D(ki,k 2 ) is a linear superposition of 
the particular solutions 

D mj (ki,k 2 ) = <k 1 ,k 2 |/,mi). (87) 

In this section we show that a general ansatz for Do(ki, k 2 ) may be obtained in terms of 21 + 1 unknown functions of 
the moduli hi, fc 2 of ki, k 2 , and of the angle 9 € [0, tt] between ki and k 2 . We also consider constraints on these 21 + 1 
functions imposed by the invariance by parity, and by the exchange (bosonic or fermionic) symmetry of D(ki,k 2 ). 



A. Ansatz for an angular momentum / 



Consider a general rotation TZ of SO(3). In the functional space in which D(ki,k 2 ) lives, it is represented by the 
unitary operator R, Rr 1 = R\ such that 

R\k u k 2 ) = \Kk u Kk 2 ). (88) 

Within the degenerate eigenenergy subspace of angular momentum I, the matrix elements of R form the unitary 
matrix p of the irreducible representation of TZ of spin /: 



{l,mi\R\l,m[) = p mim >. 



(89) 
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This matrix p is well known from group theory, where it is expressed by means of the Euler representation of TZ [40| : 

n = K z (5 1 )K y (6 2 )K z (5 3 ), (90) 

with 82 € [0, 7r], < Si, 83 < 2ir, and TZ U (4>) in SO(3) is the rotation of angle <fi around axis 2/, v = x, y, z. 
A well-known fundamental relation is then 

D(ftk 1 ,ftk 2 )=p*D(k 1 ,k 2 ) = V _1 D(ki,k 2 ) (91) 

where the spinor D is the 21 + 1 vector of components D mi , p* is the complex conjugate of the matrix p, ' p is its 
transpose, p~ 1 its inverse [45j . For fixed non-zero wave vectors kj, k 2 , one then choose the matrix 7?. such that 

ftki = fcie x (92) 
lZk.2 = k2(cos9e x + sin9e y ), (93) 

where is the non-oriented angle between ki and k 2 , 

cos# = ^T^, 9e [0,vr], (94) 

KlK 2 

and e„ is the unit vector along axis is, v = x,y, z. This choice, together with the fact that TZ is a symmetric matrix 
of determinant unity, implies 

^(k 1 Ak 2 ) = ||k 1 Ak 2 ||e z . (95) 

Remarkably, from Eq. (|9ip we see that the spinor D for arbitrary ki , k 2 is totally determined if one knows the spinor 
fW which is a function of three real variables only, 

(cos8e x +sm6e y )]. (96) 

Eq. ([9"Tj) indeed gives 

D(k 1 ,k a )=*pf«(Ai J fc 2 ,fl). (97) 

To make the ansatz more explicit, we set Z3(ki,k 2 ) = Z? m j =0 (l t i, k 2 ) without any loss of generality, so that 

D(ki, k 2 ) = Y^m =-l Pmiofmi(ki,k2,0). Using the relation (8.6-1) of we find the expression of the matrix element 
PmiO in terms of the spherical harmonics Y"/"' (9, (f>) where 9 is the polar angle and <f> the azimuthal angle. This leads 
to 

D(k 1 ,k 2 )=r^ I V /a J2 \Yr i frM*f%faM. (98) 

The last step is to calculate the angles 81 and 5 2 that enter the Euler representation of TZ. From Eq. (|90[) it is 
apparent that Si and 82 are the azimuthal and polar angles of the vector 7Ze z in the spherical coordinates of axis z. 
The Cartesian coordinates (X, Y, Z) of that vector are readily evaluated using e„ • 7Ze z — e z ■ TZ~ l e v , v = x, y, z, and 
the identities Eqs. (|92I93I95[) after their multiplication by the matrix IZ^ 1 : 

X = sin<52Cos<5i = — — - (99) 
ki 

Y = sin 8 2 sin 8 X = „ cos 9^-^] /sin 9 (100) 

V k 2 h J 

Z = cos5 2 = e z -^4^. (101) 
||ki A k 2 |] 

We note that, in the most common case I = 0, one has Y®(6, <f>) = l/\^n so that the general ansatz reduces to the 
usual form Z?(k!,k 2 ) = fn°\ki,k2,6). 
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B. Inclusion of parity and exchange symmetry 

Let us apply the parity operator to the function D (ki, ka). The change (ki, ka) ( — ki, — ka) transforms (X, Y, Z) 
into (—X, —Y, Z) so that the angle 82 is unchanged whereas 8\ is transformed into 8\ + it modulo 2ir. This changes 
Y™ 1 (82,81) into (— l) m iy ; m ' (S 2 , 81). This allows to decouple the even mi terms of Eq. (|98p. that correspond to an 
even-parity T] = +1 component of D, from the odd m; terms of Eq. (|98|) . that correspond to an odd-parity 77 = — 1 
component of D. E.g. if one considers the ansatz of I = 1 with 77 = +1, one has Y®(Q,<j)) — [3/(47r)] 1 / 2 cos 9 so that 

C(k 1 ,k 2 )=e 2 .^^/( 1) (fc 1 ,fc 2 ,0). 

Let us now consider the bosonic or fermionic symmetry requirement. Under the exchange (ki,k2) — > (ka,ki), we 
find that (X, Y, Z) is transformed into (X', Y', Z') with 

X' = X cos 9 + Y sin 9 (102) 
Y' = X sin 9 - Y cos 9 (103) 
Z' = -Z. (104) 

Then the polar 8' 2 and azimuthal <5^ angles of (X' , Y', Z') are given by S[ = 9 — Si modulo 2ir and 8' 2 = it — S 2 . From 
the following properties of the spherical harmonics, 

iHTr-M + Tr) = (-l)'lT(M) (105) 
Yf^V,® = (-l) m Yr*(9,cj)) (106) 

= Yr(9,~4), (107) 

we obtain Y^iir — 82, 9 — 81) = (— \) 1 e lme Y l ~ m (82, 8\). The bosonic (e = 1) or fermionic (e = —1) symmetry imposes 
D(k2,ki) = e£>(ki,k2) so that the unknown functions fm\ in the ansatz must satisfy 

fWfo, h,9) = e(-l) l e im ' 9 f^ mi (h,k 2 , 9). (108) 
A more convenient form is obtained by the change of unknown functions, 

fm l (h,k 2 ,6) = e-" n ' e / 2 f mi (k 1 ,k 2 ,9). (109) 

The exchange symmetry then imposes 

/*,,(*&, = e(-l) l f- mi {k 1 ,k 2 ,6), (110) 

which easily allows to restrict the calculation to the domain &2 > fci- 

By inserting the general ansatz (|9"8")l into the integral equations Eq. (I5D1) or Eq. (|84p . one may obtain general 

equations for the /m|. This was done in Ref. (3~ij in the particular case of the 3+1 fermionic problem with zero-range 
interactions, infinite scattering length and at zero energy. We shall not pursue this general issue here. 



V. CONCLUSION 



In the frame of a quite realistic two-channel model for the interaction Hamiltonian close to a magnetic Feshbach 
resonance, with three physical parameters and a potential range, we have shown that, for two relevant problems in 
the current few-body physics with cold atoms (four bosons, or three same spin state fermions plus an extra particle) 
it is possible to derive exactly simple integral equations for an unknown function of six real variables. Furthermore, 
using the rotational symmetry of the Hamiltonian and its invariance by parity, for an angular momentum Z, we have 
shown that the problem may be reduced to integral equations for I or I + 1 unknown functions of 3 real variables only, 
two wave vector moduli and one angle. 

Within the considered model Hamiltonian, this makes the exact numerical solution of the four-body problem 
tractable. We presented elsewhere a first application of this formalism in the limiting case of a zero-range interaction 
(34| . We expect that other applications will appear, in particular in the regimes where the zero-range limit with no 
extra three-body parameter is ill-defined, as it is the case for four bosons, or for the 3 + 1 fermionic problem with a 
fcrmion to extra particle mass ratio larger than ~ 13.607. 
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fc 2 abg-R*)|/|i + fc-R, | <C 1. If ibg = O(b) —> 0, the supremum is zero. If ab g is so large and negative that fc typ (— abgi?,) 1 ^ 2 > 1, 
the supremum is infinite and the condition is violated. 

As a reminder, the solution of a free space problem involves two aspects, Schrodinger's equation itself and boundary 
conditions that one must impose on the solution to recover the right physics. Implementing these boundary conditions is 
performed at two levels: (i) in the derivation of integral equations from Schrodinger's equation, when one divides by an 
energy difference E — E{ Tac , where -Ef ree is usually the unperturbed energy of some configuration, e.g. plane waves, one has 
to be careful if E — Ef Tee can be zero, in which case the resulting vanishing denominators have to be regularized by the 
inclusion of a vanishing, usually positive imaginary part, which amounts to replacing E with E + i0 + ; this division by a 
vanishing quantity also means that an additive free wave solution has to be included, that corresponds physically to an 
incoming wave, (ii) Even if all energy denominators E — E{ rcc are strictly negative, as it is the case in this paper, one in 
general still has to choose an ansatz for the unknown functions [here the functions A, B, C] that may be the sum of an 
incoming bit and of a scattered bit. E.g. the scattering of two dimers of opposite incoming wavevectors ±ko may take place 
at total negative energy, so that the integral equations derived here strictly hold [level (i) may be skipped], but still ansatz 
have to be given, C(k) = i(27r) 3 [<5(k — ko) + 5(k + ko)]p c ioscd + C SC att(k), where p c iosed is the probability for a dimer to be 
in the closed channel and C sca tt(k) is the scattered wave contribution, etc. 

One can distinguish three types of permutations. Type one includes four permutations: It corresponds to <t{1,2} = 
{3,4} (which implies <r{3, 4} = {1,2}). These four permutations give equal contributions. Type two also includes four 
permutations: It corresponds to <t{1, 2} = {1,2} (which implies <r{3, 4} = {3,4}). These four permutations give equal 
contributions. The third type includes sixteen permutations: It corresponds to the cases where there exists a unique 
i 6 {1,2} such that a(i) £ {3,4} (which implies that there exists a unique j £ {3,4} such that a(j) £ {1,2}). One can 
reduce to the case a(i) = 3, exchanging the integration variables k3 and k4 if necessary. Let i' be the element of {1,2} 
different from i. Then two cases remain to be distinguished, the case o(i') = 1 and the case er(i') = 2, which give different 
contributions. 



In particular, comparison of Eqs. (10,12) of [2J| indicates that a square is missing in the factor 7r of Eq. (10) in [24| . 

We recall that Eq. (I91|l may be derived by calculating (ki, k2|-R t \l, mi) in two different ways. First, we take to act on 
the left and we use the definitions (|88p and (JS7J), so that (ki, k2|-R t \l, mi) = D mi (IZki , 7?.k2). Second, we insert the closure 
relation 5^ m ' I'j m 'i)(^i m i\ t° the left of Ft 1 in the quantity (ki, k.2\R^ \ l,mi) , and we use Eq. (189[) . 




